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(U ; ABSTRACT 

pL^ . We present evidence that there is an M dwarf problem similar to the previously iden- 

tified G dwarf and K dwarf problems: the number of low-metallicity M dwarfs is not 
sufficient to match simple closed-box models of local Galactic chemical evolution. We 
estimated the metallicity of 4141 M dwarf stars with spectra from the Sloan Digital 
Sky Survey (SDSS) using a molecular band strength versus metallicity calibration 
developed using high resolution spectra of nearby M dwarfs. Using a sample of M 
dwarfs with measured magnitudes, parallaxes, and metallicities, we derived a relation 
that describes the absolute magnitude variation as a function of metallicity. When we 
examined the metallicity distribution of SDSS stars, after correcting for the different 
volumes sampled by the magnitude-limited survey, we found that there is an M dwarf 
problem, with the number of M dwarfs at [Fe/H] ~ —0.5 less than 1% the number 
at [Fe/H] — 0, where a simple model of Galactic chemical evolution predicts a more 
gradual drop in star numbers with decreasing metallicity. 

Key words: stars: abundances - stars: late-type - stars: statistics - Galaxy: abun- 
dances - Galaxy: evolution - Galaxy: stellar content 
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1 INTRODUCTION 

The compositions of stars provide a test for models of Galac- 
tic chemical evolution, with the general assumption being 
that the photospheric chemical composition of most stars 
represents the local Galacti c che mical composi tion w here 
they formed, van den Bergh (| 19621 ) and Schmidt l| 19631 ) first 
noted that the ratio of low-metallicity G dwarfs to solar- 
metallicity G dwarfs in the local neighbourhood is too small 
to be explained by the 'Simple model' of Galactic chemi- 
cal evolution. The Simple model, as summarized by Tins- 
ley 1 1980.1 . is a model assuming (1) the solar neighbour- 
hood can be modeled as a closed system; (2) it started as 
100% metal-free gas; (3) the initial stellar mass function 
(IMF) is constant; and (4) the gas is chemically homoge- 
neous at all times. Each of the assumptions in the Sim- 
ple model are demonstrably false, but it is an important 
starting point for the development of more complex models. 
The process of solving the G dwarf problem and explaining 
the Galaxy's chemical enrichment history has been under- 
way for nearly five decades. The solution will require dis- 
carding one or more the Simple mod el's assumptio n s. The 
intro d uction of a variable IMF (e., 
1 19961 : iMartinelli fc Matteuccil l200i 



ISchmidt 1963 



iRomano et al 



variable star forma t ion rates with intermittent mixing (e.g. 
iMalinie et "al]|l993l : [Caimmill2008l ). and / or inflow or out- 
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flow of material fe.g. lWvse fc Gilmorelll995l : |Pagelll200ll ) into 
models of Galactic chemical evolution produces metallicity 
distributions in better agreement with observations. 

The G dwarf problem does not apply only to the local 
neighbo urhood: the G dwarf problem exist s in other galaxies 
as well (|Worthev. Dorman. fc Jon"eill996l ). While G dwarf 
lifetimes would suggest that some G dwarfs may have left 
the main sequence during the Galaxy's lifetime, the paucity 
of low-metallicit y stars locally also extend s to the longer- 
lived K dwarfs l|Casuso fc BeckmanI liooi ) . Mould (Il978l ) 
suggested that there may be an M dwarf problem: using 
spectra of six M dwarfs and infrared spectroscopy of sixteen 
old disc M dwarfs, he estimated the abundance dispersion 
for M dwarfs and found evidence that suggested a common 
chemical history of G and M dwarfs. 

We have directly tested the existence of an M dwarf 
problem by estimating the metallicity of 4141 M dwarf stars 
with temperatures in the range where our analysis method 
is valid, 3500K T^s ^ 4000K, using CaH and TiO molec- 
ular band strengths measured from Sloan Digital Sky Sur- 
vey (SDSS) spectra (|York et al.. .2000) and find that there 
is an M dwarf problem in the Galactic disc similar to the 
previously known G and K dwarf problems. This does not 
suggest a new solution to the M, K, and G dwarf problems, 
but rather shows that for all stars with main sequence life- 
times comparable or larger than the age of the Galaxy, and 
where chemical abundance surveys have been performed, the 



© 0000 RAS 



2 Vincent M. Woolf and Andrew A. West 



number of low-metallicity stars is insufficient to matcfi the 
Simple model. 



2 ESTIMATING METALLICITIES 

M dwarfs are inherently faint objects, and therefore must be 
nearby to appear bright enough to obtain high-resolution 
spectra, A/AA > 30,000, of sufficient quality to measure 
atomic absorption lines for abundance analyses, while using 
a reasonable amount of telescope time. 

Molecular band indice s CaH2, CaH3, and Ti05 
l|Reid, Hawlev. fc GizidlToOSI ) measured from low-resolution 
spectra, A/AA ~ 1800, have been shown to be useful in clas- 
sifying cool dwarfs. CaH2 and CaHS correlate well with spec- 
tral type and the combination of CaH2 or CalfS with Ti05 
separates cool dwarfs into rough m etallicity classes: dwarfs. 



subdwarfs, and extreme subdwarfs |Gizis| 19971 : iLepine et al.l 
I2OO3I ). Woolf, Lepine, & Wallerstein I2009I ) used metallicities 
of 88 M dwarfs, estimat ed using atomic lines measured f rom 
high-resolution spectra (| Woolf fc Wallerstein|[2005l . l2006l ). to 
develop and calibrate a method to use CaH2, CaHS, and 
Ti05 molecular band indices to estimate M dwarf metal- 
licities more precisely. CaH2, CaHS, and their sum (CaH2 
-I- CaHS) were found to correlate well with effective tem- 
perature (Toff) a nd Ti05 was found to vary with both Tcff 
and metallicity l|Woolf fc WallersteinI I2OO6 1 The three in- 
dices were used with the calculated [Fe/H j metallicites of 
the 88 M dwarfs to derive the relation 



[Fe/H] ^a + biCT.o/CaH), 



where a 



-1.685 ± 0.079, b 



(1) 

1.6S2 ± 0.096, and the 
metallicity index Cxio / CaH is defined as described in Lepine, 
Rich, fc Shara (|20o3). This calibration of the molecular 
band index versus metallicity relation is valid for stars with 
S500K s; Teff 4000K and -1.5 < [Fe/H] < +0.05. As de- 
scribed in lWoolf et all l|2009l '). the [Fe/H] versus Cxio/CaH fit, 
should be accurate to ±0.S dex for stars with [Fe/H] > —1.0 
and to ±0.5 dex for stars with [Fe/H] < —1.0. based on 
the maximum errors the empirical data appear to allow, 
as opposed to standard statistical errors based on repeated 
measurements of the same quantity. Most errors should be 
less than half of those maxima. 

We used CaH2, CaHS, and Ti05 indices measured from 
SDSS spectra to estimate metallicities for the 4141 M dwarfs 
in our sample. The stars were placed in 0.1 dex [Fe/H] bins 
from —1.5 to 0.1. This extrapolates slightly beyond the lim- 
its for which our method was calibrated. The number of stars 
in bins with [Fe/H] less than —0.50 is very small, so extrap- 
olating to [Fe/H] = —1.55 at the low metallicity end makes 
little difference. The number of stars in the [Fe/H] ~ -fO.l 
bin is less certain because the entire bin is outside our cal- 
ibration, but the bin is included to show that, if our cali- 
bration continues to be correct a little beyond where testing 
stopped because of a lack of stars for the calculation, it ap- 
pears that the trend in number of stars versus metallicity 
drops for [Fe/H] > 0. 



wliere [X] = Iogio(X)s 



logio(X)sun- 



3 STELLAR DATA SET 

Our spectroscopic sampl e contains 4141 lo w-mass stars se- 
lected from the SDSS (|York et all I2OOOI ') Munn 'special 
plates,' which were made public as part o f the SDSS 
Data Release 4 jAdelman-McCarthv et aLlbOOd V The Munn 
plates were targeted to be a magnitude-limited sample of 
M and late-K dwarfs with i < 18.26, i — z > 0.2 and be 
located near the / = 12S°, b — — 6S° Galactic sight line. 
This sample was previously used to study the kinematics of 
the lo cal Milky Way thin and thick discs (iBochanski et al] 
|-JiMi7;J i to roughly 2 kpc above the Galactic plane and are 
contained in the larger SDSS M dwarf spectroscopic samples 
jWest et al. 2004, ,2008, 10 U). We analysed the 8859 spec- 
tra that were observed us ing the Hammer spectral typing 
facility (ICovev et al.ll2007l ) and rejected stars that had poor 
quality (SNR < S near Ha) and spectral types earlier than 
K5 (resulting in 7714 stars). 

As part of our analysis, we computed radial velocities 
(RV ) by cross-correlating ea ch spectrum with the appropri- 
ate IjBochanski et al. I l2007bl ) M dwarf template. In the case 
of K5 and K7 dwarfs, we used the MO template for our RV 
determinations. The computed velocities were used to cor- 
rect the stellar spectra to zero radial velocity. We then com- 
puted the Ti05, CaH2 and CaHS molecular band indices 
(|Reid et all Il995l : ICizid Il997l ) from the velocity-corrected 
spectra. Additional stars were removed from the sample if 
their molecular band indices, temperature, or [Fe /H] fell out- 
side the range covered by the Woolf et al. (2009?) metallicity 
calibration, leaving the 4141 stars for our study. All stars 
with temperature and [Fe/H] within the metallicity calibra- 
tion have i — z colours larger the 0.2 value used in the initial 
selection, so the colour limit should introduce no selection 
effects for our metallicity statistics. 



4 VOLUME AND STELLAR NUMBER 
DENSITY CORRECTIONS 

4.1 Volume correction for V and R magnitudes 

Low metallicity main sequence stars (subdwarfs) are less lu- 
minous than higher metallicity sta rs with the same temper- 
ature or spectral tvpe (jRciz. 19541: iMould fc McElrovl Il978l : 
lAke fc Greensteiiilll980l : Bochanski et al.ll2011 ) . Because the 
SDSS M dwarf sample was selected to fall within a magni- 
tude range, the subdwarf M stars must be closer on average 
than the higher metallicity stars, and thus sample a smaller 
volume of our Galactic neighbourhood. We must correct for 
this to avoid introducing a Malmquist Bias-like effect and 
overcounting the numb er of more luminous stars. 

We used stars from I Woolf et a l. (2009) for which paral- 
lax data were available and thus absolute magnitudes could 
be calculated to find the luminosity variation with metal- 
licity for M dwarfs in our temperature range. An absolute 
magnitude difference function was found by fitting a three 
dimensional 'surface' to the [Fe/H], CaH2-|-CaHS, and My 
data points calculated for the stars. CaH2-|-CaHS acts as a 
temperature proxy in the calculated fit. We found that for 
the best-fitting surface, the magnitude change with metal- 
licity did not depend on CaH2-f CaHS (temperature) within 
the range of our metallicity calibration. For example, two 
S500 K stars with [Fc/H] = 0.0 and -1.0 differ in magnitude 
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Figure 1. Absolute V magnitudes for early M dwarfs calculated 
using Equation [2] versus those calculated using observed V mag- 
nitudes and parallax measurements. The diagonal line indicates 
where points would fall if the model and observed values were 
identical. 



by the same amount as two 4000 K stars with [Fe/H] — 0.0 
and -1.0. The fit to the data points is 

Mv = a + b[Fe/H] + c[Fe/H]^ + d(CaH2 + CaH3) (2) 

+e(CaH2 + CaH3)^ 

where a = 25.356, b = -0.6394, c = 0.8455, d = -15.778, 
and e — 3.250. The root mean squared deviation between 
the absolute V magnitudes calculated using parallax and ap- 
parent V magnitudes and those calculated using the above 
formula is 0.450. The My values calculated from Equation 
[21 Mv (model), are compared to those calculated from ob- 
served V magnitudes and trignometric paralax, Mv (ob- 
served), in Figure[T] We did a similar fit using the absolute R 
magnitudes and found coefficients a = 24.334, b — —0.8258, 
c = 0.6964, d = —16.379, and e = 3.696, with the root mean 
squared deviation between the Mr values calculated from 
paralax and R measurements and from the formula being 
0.485. 

Because there is no temperature dependence in the dif- 
ference in magnitude (or luminosity) for M dwarfs of the 
same temperature but with different metallicities there is 
also no temperature dependence in the volume correction 
due to metallicity. We can apply the same volume sampling 
correction for all stars in our sample without considering 
temperature. 

The magnitude-limited volume sampled for each metal- 
licity bin is proportional to the distance cubed. The distance 
to a star is proportional to the square root of its luminosity 
for a given apparent magnitude. So for a given stellar lumi- 
nosity, L, the sampled volume is proportional to L^^^. Sub- 
dwarfs with smaUer luminosities are sampled from a smaller 
volume, so we multiply the numbers of stars in the metallic- 
ity bins by the inverse of the volume factor, L"^''^, to correct 
for the volume difference. The V and R luminosity and abso- 
lute magnitude differences and the volume correction factor 
for di fferent [Fe/H] values as calculated using Woolf et al.l 
l|2009h stars are shown in Table [T] 



4.2 Volume correction for SDSS i magnitude 

The SDSS stars in our sample were selected for observa- 
tion using an SDSS i magnitude cutoff. We must therefore 
use AMi, the difference in i absolute magnitude caused by 
varying metallicity, in our volume correction. 

None of the stars used for our volume correction due to 
metallicity differences have SDSS photometry available, so 
they cannot be used to calculate AA/i directly. Instead, we 
first use AMv and AM_r to estimate AMr, the difference in 
SDSS r magnitude due to metallicity differences. We then 
use A(r — i), the difference in r — i colour due to metal- 
licity differences, to calculate AMi to use for our volume 
correction. 

The SDSS r photometric band overlaps the V and R 
bands. Differences in metallicity should therefore cause dif- 
ferences in r absolute magnitude similar to the differences 
found for V and R. The corrections for V and R differ at 
most by 0.06 magnitudes. Assuming that AMr is equal to 
the mean of AAfv and AM_r should therefore introduce er- 
rors smaller than a few percent. 

We have therefore estimated the necessary correc- 
tions using updated synthetic spectra (Hauschildt, pri- 
vate communication) based on PHOENIX calculations 
(Hauschildt, Allard, fc Baron 1999 ), finding the mean dif- 
ferences in r — i colours in the spectra relative to the [Fe/H] 
= 0.0 spectrum for temperatures between 3500 K and 4000 
K. We then calculated the difference in i magnitude due to 
metallicity differences: AMi = AM,- — A(r — i). The cor- 
responding i luminosity ratio and volume corrections were 
calculated in the manner described previously for V and R 
and are reported in Table [T] 

An alternate method to calculate A(r — i) would use a 
published empirical colour-temperature calibration to estim- 
age B — V, V — R, and V — I colours for stars of M dwarf tem- 
peratures then use an existing empi rically-derived f o rmula 
to transform these to an r — i colour. IWorthev fc Led ()201lh 
provide such a colour-temperature calibration. The M dwarf 
metallicities used in their calibration are based on kinemat- 
ics, with young disc stars assigned [Fe/H] = —0.1, old disc 
stars —0.5, and halo stars —1.5, unless the star is in a clus- 
ter, in which case the cluster metallicity was used. If we use 
their empirical calibration rather than the colours from syn- 
thetic spectra, we find much smaller corrections for A(r — i). 
For example, the r — i correction for [Fe/H] = —1.5 stars is 
—0.3 if calculated from synthetic spectra, but is about 0.02 
if calculated from the empirical calibration data. 

Our goal is to compare the trend for M dwarf stars with 
a model which has predicted more low-metallicity stars than 
have been observed. We will therefore avoid the possibility of 
under-counting low-metallicity stars by using the larger cor- 
rection based on synthetic spectra which tends to increase 
the corrected number of low-metallicity stars. In the end, 
because the number of observed stars drops off so quickly 
with decreasing metallicity, our conclusions do not depend 
on which of these corrections we use. We correct by multi- 
plying the number of stars in each metallicity bin by the val- 
ues in the volume"^ (i) column in Table [1] The empirically- 
based difference between AMv and AMi is very small, so if 
someone prefers not to use the correction based on synthetic 
spectra, one can multiply by the value in the volume"^ (V) 
column instead. 
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We note that the SDSS sample has both a faint i mag- 
nitude hmit and a bright i magnitude hmit, meaning that 
the volumes sampled are effectivly spherical shells. The ra- 
tio of the distances to the faint magnitude cutoff is the same 
as the ratio of the distances to the bright magnitude cutoff 
however, so the volume ratio of the unsampled inner spheres 
is the same as the volume ratio of the spheres defined by the 
faint limit: no additional correction is required to account 
for the bright magnitude cutoff. 



4.3 Correction for stellar number density 
differences 

The observational direction for our SDSS M dwarf sample, 
near I = 123°, b — —63°, means that the more distant stars 
in the sample are farther from the plane of the Galactic disc 
and thus in a region with a lower stellar number density. 
Without correcting for this effect we would overestimate the 
fraction of fainter, low-metallicity subdwarf stars which are 
found in a denser region of the disc. We use the difference in 
distances due to metallicity and a model of the Milky Way 
stellar numb er densit y distribution to correct for this effect. 

.Bochan ski et al. I lloii) report absolute r magnitudes, 
Mr, for M dwarfs based on photometric parallax. Figure 2 
from their paper plots Mr vs spectral type. From this fig- 
ure we find that the typical Mr for stars in our data set, 
mostly MO to M2 dwarfs, is about 9.0. This is consistent 
with the absolute magnitude of the few warmer M dwarfs 
with reported SDSS ph otometry and trigonometric parallax 
l|Davenport et aLlbOOd ). A Mr = 9.0 star with r = 17.5, the 
mean observed r magnitude for our sample, is at a distance 
d = 501 PC. 

lJuric et al.l (|2008l ) report the stellar number density dis- 
tribution for the Milky Way, providing an exponential equa- 
tion to calculate estimates of changes in density with radial 
distance from Galactic centre and perpendicular distance 
from the plane of the disc (Juric et al. Equation 23). Stellar 
number density variation due to differences in radial dis- 
tance from Galactic centre is minimal for our stars, given 
the direction and distance to the sample. We assume the 
radial terms from the Juric et al. equation are constant and 
recognize that Zq is above the plane of the Galaxy and the 
direction to the stars in our sample places them below the 
plane. The equation thus becomes p = C exp(— — ), where 
C is the stellar number density at the vertical centre of the 
disc at the Solar radial distance, Zq is the Solar distance 
above the Galactic plane, Z is the distance to the star mea- 
sured in the direction perpendicular below the plane, and 
H is the scale height for the disc, thin or thick. Note that 
this form of the equation is not valid for stars above the 
plane or closer to the centre of the plane than the Sun. We 
use this modified equation to find the ratio of the stellar 
number density p for a star in our sample with a distance Z 
compared to the local density : 

— =exp( — ). 



(3) 



lJuric et all (|2008l 'l report H = 300 pc for the thin disc, H = 
900 pc for the thick disc, and Zq = 20 pc. 

To make a rough correction for differences in stellar 
number density caused by metallicity we first use AMr from 
Table [T] to find the luminosity ratio relative to a [Fe/H] = 



[Fe/H] 

Figure 2. Fraction of stars found in each 0.1 dex [Fe/H] bin. The 
daslicd line represents tlie raw data. The solid line represents the 
volume and stellar number density corrected data. The dotted 
curve is the distribution predicted by the Simple Model. 

star for each metallicity bin; I/Iq = 100"'^*^''/''. The ratio 
of the distances for stars with the same apparent magnitude 
is then cL/cLq = ^JTfk). If we then assume that the mean 
distance to a [Fe/H] = star in our sample is 501 pc we can 
use the distance ratio to find the mean distance to stars in 
the other metallicity bins. Because h ~ —63° for our sample, 
the distance perpendicular to the Galactic plane vs, Z = 
d sin 61° - Zq. 

We calculate the ratio of the stellar number density 
at the typical perpendicular distance, Z, for stars in each 
metallicity bin relative to the local Solar neighborhood den- 
sity using Equation [3l We calculate the ratio for both thin 
disc stars, (^)thin, and thick disc stars, (^)thick, and use 
these and the local 12% thick disc fraction l|juric et al.l 
|2008| ) to calculate the density ratio ^ = 0.88(^)thin + 
0.12(^) thick- For each metallicity bin we divide this den- 
sity ratio by the value found for the [Fe/H] = bin 
to get where po is the stellar number density at the 
typical distance to solar metallicity stars in our sample. To 
correct for differences in stellar number density we multiply 
the number of stars in each metallicity bin by the inverse of 
this density ratio, as reported in Table [T] We note 

that the volume correction and the stellar number density 
correction act in opposite directions; fainter, lower metallic- 
ity stars are sampled in a smaller volume, but in a region 
with a higher density of stars. 



5 M DWARF METALLICITY TREND 

The raw and corrected fractions of M dwarfs in our sample 
in each metallicity bin are shown in Table [2] and in Figure 
(2] The distribution is centred at about [Fe/H] = 0.0 and has 
a gaussian full width half maximum (FWHM) of 0.30 dex 
and a standard deviation (dispersion) of 0.13 dex. 

Other studies have found that the number of 
stars peaks at a metallicity centred in the range 
0.25 < [Fe/H ] < 0.0 for G and K dwar fs (e.g . 



Wvse fc Gilmord 1 19951; iRocha-Pinto fc MacidI 



Favata. Micela. fc Sciortiiiol 1997 : J0rgensen 



1996 : 



2000; 
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Table 1. M dwarf absolute magnitude differences and luminosity, volume, and stellar number density cor- 
rection factors for sampling differences compared to Solar [Fe/H]. Values in the volume" -"^ (i) and {'^)~^ 
columns are the factors by which the numbers of stars in each [Fe/H] bin should be multiplied to correct for 
volume sampling and stellar number density differences. A(r — i) is the colour difference relative to [Fe/H] 
= 0.0 due to metallicity difference. 



fFc /HI 


AMv 


AMr, 


L(V) 




volume ^ {V) 




Afr — i) 


AAfi 


Ui) 


volume ("i) 




+0.1 


-0.056 


-0.075 


1 


056 


1 


072 


0.925 


-0.066 


+0.050 


-0.116 


1.112 


0.852 


1.038 


0.0 


0.000 


0.000 


1 


000 


1 


000 


1.00 


0.00 


0.000 


0.00 


1.00 


1.00 


1.00 


-0.1 


0.072 


0.090 





936 





920 


1.10 


0.081 


-0.037 


0.118 


0.897 


1.18 


0.957 


-0.2 


0.161 


0.193 





862 





837 


1.25 


0.177 


-0.074 


0.251 


0.794 


1.41 


0.909 


-0.3 


0.268 


0.311 





781 





751 


1.45 


0.289 


-0.111 


0.400 


0.692 


1.74 


0.859 


-0.4 


0.391 


0.442 





698 





656 


1.72 


0.416 


-0.148 


0.564 


0.595 


2.18 


0.807 


-0.5 


0.531 


0.587 





613 





582 


2.08 


0.559 


-0.184 


0.743 


0.504 


2.79 


0.756 


-0.6 


0.688 


0.746 





531 





503 


2.59 


0.717 


-0.203 


0.920 


0.429 


3.56 


0.707 


-0.7 


0.862 


0.919 





452 





429 


3.29 


0.890 


-0.222 


1.112 


0.359 


4.65 


0.659 


-0.8 


1.052 


1.105 





380 





361 


4.28 


1.079 


-0.240 


1.319 


0.297 


6.18 


0.615 


-0.9 


1.260 


1.306 





313 





300 


5.70 


1.283 


-0.258 


1.541 


0.242 


8.41 


0.573 


-1.0 


1.485 


1.521 





255 





246 


7.78 


1.503 


-0.277 


1.780 


0.194 


11.7 


0.535 


-1.1 


1.726 


1.749 





204 





200 


10.9 


1.738 


-0.284 


2.022 


0.155 


16.3 


0.501 


-1.2 


1.985 


1.991 





161 





160 


15.5 


1.988 


-0.291 


2.279 


0.123 


23.3 


0.470 


-1.3 


2.260 


2.247 





125 





126 


22.7 


2.253 


-0.299 


2.552 


0.953 


34.0 


0.442 


-1.4 


2.552 


2.518 


0.0953 


0.0984 


34.0 


2.535 


-0.306 


2.841 


0.0730 


50.7 


0.418 


-1.5 


2.861 


2.802 


0.0717 


0.0757 


52.1 


2.831 


-0.313 


3.144 


0.0552 


77.0 


0.397 



Table 2. Numbers and sample fraction of M dwarfs at given metallicities. The last two columns are cor- 
rected for different volumes, volume~^(i), and stellar number densities, being sampled for different 
metallicities. 



[Fe/H] 


number 


fraction 


corrected number 


corrected fraction 


+0.1 


1318 


0.3183 


1166.3 


0.2685 


0.0 


1491 


0.3601 


1491.0 


0.3433 


-0.1 


868 


0.2096 


977.7 


0.2251 


-0.2 


312 


0.0753 


401.0 


0.0923 


-0.3 


99 


0.0239 


147.9 


0.0340 


-0.4 


33 


0.0080 


58.1 


0.0134 


-0.5 


8 


0.0019 


16.9 


0.0039 


-0.6 


3 


0.0007 


7.75 


0.0018 


-0.7 


1 


0.0002 


3.07 


0.0007 


-0.8 


2 


0.0005 


7.60 


0.0018 


-0.9 


1 


0.0002 


4.82 


0.0011 


-1.0 


3 


0.0007 


18.8 


0.0043 


-1.1 














-1.2 














-1.3 














-1.4 


2 


0.0005 


42.3 


0.0097 


-1.5 















2004 



2004 



Rocha-Pinto et all I2OO0I: lAUende Prieto et all 
Casuso fc BeckmanI l2004l : | Nordstrom et al.l 
Luck fc Heiteil 120051 ). We find that the peak for M 
dwarfs is at the high end of this range and note that our 
peak is more narrow than that typically found. There 
has been some "narrowing" of the reported metallicity 
distribution over time since the first reports of the G dwarf 
problem in the 1960's, probably as a result of a reduction of 
the uncertainties as the methods for calculating elemental 
abundances have improved. 

Our results for M dwarf s most closely match those 
found by Favata et al. (|l997l ') for K dwarfs, however. In 
their study of G and K dwarfs, they separated stars into 
two groups with effective temperatures hotter and cooler 



than 5100 K and found very different metallicity distribu- 
tions. They did not report mean and standard deviations 
for their corrected distributions, but we estimate from the 
values on their plots that the warmer stars showed a dis- 
tribution similar to those previously reported for G dwarfs, 
centred at [Fe/H] « -0.28 with a fa 0.29, and that the 
cooler star distribution dropped off more sharply toward 
lower metallicities and was more sharply peaked with acen- 
tre at [Fe/H] ^ -0.03 and a « 0.17. Luck fc Heiter (J200^, 
report elemental abundances for stars within 15 pc. Their 
sample contains mostly G and K dwarfs. When we sepa- 
rate their stars by temperature we find that the metallic- 
ity distribution for the cool stars, Tefi ^ 5100K, has mean 
[Fe/H] = 0.0 and a = 0.17, while the warmer stars have 
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[Fe/H] 



Figure 3. Theoretical (dashed curve) and observed (sohd line) 
cumulative stellar metallicity distributions. The theoretical val- 
ues are calculated using the Simple model. S/Si represents the 
fraction of stars that have metallicites less than Z. 



mean [Fe/H] = —0.11 and a = 0.28, again showing an identi- 
fiable diflerence with temperature. It appears that the trend 
of higher mean metallicity and narrower peaks seen for the 
metallicity distributions for cooler stars seen in the Favata 
et al. and Luck & Heiter data applies to the even cooler M 
dwarfs in our sample. 

'Simple' models of Galactic chemical evolution produce 
higher numbers of long-lived low metallicity stars than we 
find for M dwarfs and othe rs have found for G and K dwarfs. 
lAudouze fc TinsievI |l97^) use the equation 



S_ 

Si 



Z/Zi 

1^1 



Ml 



(4) 



to calculate 5*/ Si , the present day cumulative stellar metal- 
licity distribution, or the fraction of stars that have metal- 
liciies less than Z. In equation^ Zi is the present day metal- 
licity and fj,i is the fraction of the local baryonic matter that 
is now interstellar matter (as opposed to that contained in 
stars and stellar remnants). The smooth dotted curve in Fig- 
ure [5] is derived from the cumulative distribution described 
by equation U given that the increase in the cumulative dis- 
tribution must be caused by stars in each successive metal- 
licity bin. 

In Figure [3] we compare S/Si to the value cal- 
culated from our M dwarf metallicities. We calculated 
S/Si from equation 3] w here we assume that fii = 0.27 
l|Holmberg fc Flvnn|[2000l ') and let log{Zi/ZQ) = +0.1. The 
theoretical S/Si calculated for the Simple model will be 
somewhat different for other assumed Zi and pi , and a fairly 
large range of values for fii are possible, given the 5 percent 
uncertainty reported bv lHolmberg fc Flvnnl (|200G| ). But for 
all reasonable values the theoretical curve drops off much 
more slowly toward low metallicity than the observed curve 
does. 



of M dwarfs peaks at [Fe/H] = 0.0 and drops off quickly at 
lower metallicities. 

We must be careful in comparing our M dwarf num- 
bers to models of Galactic chemical evolution because our 
M dwarfs are not necessarily representative of the solar 
neighbourhood. Galactic chemical evolution models nor- 
mally include a region about 1 kpc wide in the Galactic 
plane and about 1 kpc perpendicular to the plane to in- 
clude stars that formed in the thin disc but which now leave 
the local neighbourhood l|Tinslevl .l980). The SDSS magni- 
tude lirnits_jnBafl_thB;t_^^ out to about 2 
kpc (|Adelman-McCarthv et al]|2006l ) and that no nearby M 
dwarfs are included. The part of the survey which includes 
our M dwarf spectra is centred near the southern Galactic 
cap (Bochanski ct al. 2007a). Much of the sampled region is 
more distant than the thin disc scale height, so we expect 
a larger fraction of thick disc stars to be included in the 
sample than would be found locally. 

Including a larger fraction of stars outside the thin disc 
would presumably mean including more low-metallicity stars 
than would be found in a sample from the local neighbour- 
hood. The fact that even with the probable inclusion of more 
thick disc stars we still find a paucity of low metallicity stars 
strengthens the case for the existence of an M dwarf prob- 
lem. 

In addition to our sample being made up of non-local 
stars, the fraction of thin disc versus thick disc stars should 
vary with metallicity because metallicity is correlated with 
luminosity and distance for late main sequence stars. We 
corrected for differences in sampling volume caused by low- 
metallicity stars being less luminous and thus being observed 
at smaller distances for similar apparent magnitudes. We did 
not correct for the difference in population sampling caused 
by the variation of brightness with metallicity. Nor did we 
correct for differences in average stellar metallicity with dis- 
tance perpendicular to the disc. Solar metallicity stars in our 
sample are at a larger average distance from the Galactic 
disc than low-metallicity stars with the same temperature. 
We therefore expect the solar metallicity stars to be found 
in a region with a larger fraction of th ick di sc stars and a 
smaller average metallicity. Ivezic et al. (|2008h found no cor- 
relation between metallicity and kinematics in the Galactic 
disc, but did measure a 0.2 dex variation in metallicity from 
500 pc to several kpc above the disc. Because the metallicity 
variation with distance in our sam ple should be small, less 
than the 0.2 dex variation found bv llvezic et ai] ||2008D . our 
lack of a correction for metallicity or population differences 
in our sample should not alter our results appreciably. 

We find evidence that the small fraction of low- 
metallicity stars seen for G dwarfs is observed for M dwarfs 
as well, providing another strong constraint on models of 
Galactic chemical evolution. The solution to the G, K, and 
M dwarf problem will require using models which abandon 
the assumptions of the 'Simple model' as discussed in the 
introduction. 



6 CONCLUSIONS 

Our data show that there is an M dwarf problem similar to 
the previously known G and K dwarf problems. The number 
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